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FIELD OF THE INVENTION 

The invention relates to securing of information utilizing 
imaging technologies and more specifically to phase and 
amplitude encryption and decryption of images. 

5 BACKGROUND OF THE INVENTION 

It is well known to form an image by phase contrast imaging 
methods in which phase modulation of light is converted into 
intensity modulation. As opposed to intensity modulation, 
phase modulation does not involve loss of energy. 

10 In published patent application no. WO 96/34307 which is 
hereby incorporated by reference, a phase contrast imaging 
method is disclosed for calculating phasor values of a phase 
mask for synthesizing a desired intensity pattern. 

In "Securing Information with Optical Technologies", Bahram 
15 Javidi, Physics Today Vol. 50, No. 3, March 1997, pp 27-32, a 
method and system for optically securing information is 
proposed. An image is encrypted in a 4f-lens optical 
configuration (i.e. comprising two Fourier transforming 
lenses) by Fourier transforming the image and an input phase 
20 mask with a first lens of the 4f-lens system. A Fourier plane 
phase mask having phasor values e ib(a 'P } is positioned in the 
Fourier plane of the first lens and an encrypted image is 
formed with a second lens that Fourier transforms images in 
the Fourier plane of the first lens. The encrypted image is 
25 decrypted in a similar 4f-lens configuration in which another 
Fourier plane phase mask (the key) having phasor values 
e -ib(a,P) ig positioned in the Fourier plane of the first 
lens. 

It is a disadvantage of the known cryptographic method that 
30 encrypting an image both in the object plane and the Fourier 
plane leads to generation of speckle patterns in the 
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decrypted image thereby corrupting data having been 
encrypted. 

It is another disadvantage of the known cryptographic method 
that extremely accurate three dimensional positioning of the 
5 phase mask in the Fourier plane is required for successful 
encryption and decryption. 

It is yet another disadvantage of the known cryptographic 
method that both amplitude and phase have to be recorded in 
the encrypted mask. 

10 SUMMARY OF THE INVENTION 

It is an object of the present invention to provide an appar- 
atus of the above kind which apparatus is robust, compact, 
simple to design and relatively cheap to manufacture. 

It is another object of the present invention to provide an 
15 improved method and apparatus without the above-mentioned 
disadvantages . 

It is still another object of the present invention to pro- 
vide an improved method and apparatus for optically securing 
data utilizing phase contrast imaging. 

20 According to the invention a method is provided of decryption 
of an encrypted image having a non- encrypted image intensity 
pattern I(x',y'). 

The image is encrypted into a mask having a plurality of mask 
resolution elements (x^y,,,) by encoding the image into the 

25 mask with an encoded phase value ^(x m ,y m ) and an encoded 
amplitude value a(x m ,y m ), and by encrypting the mask by 
addition of an encrypted phase value ^(^'Vm) to tlie encoded 
phase values ^fx^y^ and by multiplication of an encrypted 
amplitude value a c (x m ,y m ) with the encoded amplitude value 

30 Mx^y^. Thus, each mask resolution element (x^y^ 
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modulates the phase and the amplitude of electromagnetic 
radiation incident upon it with the complex value 
a ix*. y m ) a c (x^, y m ) e 1 * (jm -y»> +i * c . 



The method comprises the steps of radiating electromagnetic 
5 radiation towards the mask, and inserting into the path of 
the electromagnetic radiation a complex spatial 
electromagnetic radiation modulator comprising modulator 
resolution elements (x d ,y d ) f each modulator resolution 
element (x d ,y d ) modulating the phase and the amplitude of 

10 electromagnetic radiation incident upon it with a 
predetermined complex value a d (x d , y d ) e^^dtxd^yd) ^ the 
decrypting phase value 0 d (x d ,y d ) and the decrypting amplitude 
value a d (x d ,y d ), respectively, of a modulator resolution 
element (x d ,y d ) being substantially equal to -0 C (^'Yin) and 

15 a^Mx^yro), respectively, of a corresponding mask resolution 
element (x m ,y m ) , and imaging the mask and the electromagnetic 
radiation modulator onto the image having the image intensity 
pattern I(x',y'). 



The complex spatial electromagnetic radiation modulator may 

20 be positioned anywhere in the path of the electromagnetic 

radiation. When the complex modulator is positioned adjacent 
to the mask, the encrypted mask is decrypted when the 
decrypting phase value <£ d (x d ,y d ) and the decrypting amplitude 
value a d (x d ,y d ) , respectively, of a modulator resolution 

25 element (x d ,y d ) being substantially equal to -<£ c (Xm,y m ) and 
a c ' x (x^ y m ) , respectively, of a corresponding adjacent and 
aligned mask resolution element (x^y^ . The complex 
modulator may be physically positioned remotely from the mask 
and imaged onto the position of the mask whereby the complex 

3 0 modulator is virtually positioned adjacent to the mask. The 
complex modulator may be positioned virtually or physically 
in another position than adjacent to the mask still 
performing the decrypting function when the amplitude and 
phase values of the complex modulator are transformed 

3 5 accordingly, e.g. by a Fresnel propagator. 
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It is an important advantage of the present invention that 
the encrypted image may selectively require recording of 
phase values or amplitude values or a combination thereof . 

It is another advantage of the present invention that 
5 decryption is performed in a plane adjacent to the mask or an 
equivalent plane whereby generation of speckles in the 
decrypted image is suppressed. 

It is yet another advantage of the present invention that 
there is no requirement of positioning a radiation modulator 
10 in the Fourier plane whereby an accurate three dimensional 
positioning requirement is avoided. 

It is still another advantage of the present invention that 
only a single key for encryption is required. 

It is preferred to encode the non- encrypted image into a 
15 phase mask and encrypt the phase mask by adding encrypting 
phase values to the encoded phase values. An encrypted phase 
mask is extremely difficult - if not impossible - to 
replicate by counterfeiters. Further, phase masks may be 
readily produced while masks requiring recording of both 
20 amplitude and phase are extremely complicated to produce, 

typically requiring production of two masks to be accurately 
superpositioned. The same applies to the decrypting complex 
modulator. 

The imaging may be performed with a common path 
25 interferometer, such as a phase contrast imaging system, a 

dark field imaging system, a field absorption imaging system, 
a point diffraction imaging system, a Smartt interferometer, 
a Schlieren interferometer, etc, and any combination hereof. 

According to a preferred embodiment of the invention a method 
30 is provided for synthesizing an intensity pattern with low 

loss of electromagnetic energy, comprising spatial modulation 
of electromagnetic radiation with a spatial phase mask for 
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modulation of the phase of the incident electromagnetic 
radiation by phasor values of individual resolution elements 
of the spatial phase mask, each phasor value being determined 
in such a way that 

5 1) the values of the Fourier transformed phasors attains 

predetermined values for predetermined spatial fre- 
quencies, and 

2) the phasor value of a specific resolution element of 
the spatial phase mask corresponds to a distinct 
10 intensity level of the image of the resolution 

element in the intensity pattern, 

and a spatial phase filter for phase shifting of a part of 
the electromagnetic radiation, in combination with an imaging 
system for generation of the intensity pattern by interfer- 
15 ence in the image plane of the imaging system between the 
part of the electromagnetic radiation that has been phase 
shifted by the phase filter and the remaining part of the 
electromagnetic radiation. 

Although, the present method is related to encoding of spa- 
20 tial phase masks in two spatial dimensions (planar encoding) , 
the principles of the method may be utilized for phase encod- 
ing in one to three spatial dimensions and/or in the temporal 
dimension. 

The electromagnetic radiation may be of any frequency range 
25 of the electromagnetic spectrum, i.e. the gamma frequency 

range, the ultraviolet range, the visible range, the infrared 
range, the far infrared range, the X-ray range, the microwave 
range, the HF (high frequency) range, etc. The present method 
is also applicable to particle radiation, such as electron 
30 radiation, neutron radiation, etc. 

Preferably, the electromagnetic radiation is monochromatic or 
quasi -monochromatic so that the energy of the electromagnetic 
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radiation is concentrated in a narrow frequency bandwidth. As 
the intensity pattern is synthesized by interference of two 
electromagnetic waves emitted from a common source of elec- 
tromagnetic radiation but the phases of which have been 
5 changed differently, it is required that the frequency range 
of the emitted electromagnetic radiation is sufficiently 
narrow to ensure that the two waves of electromagnetic radi- 
ation are coherent so that their superposition generates the 
desired intensity pattern. If the frequency range is too 

10 broad, the two waves will be incoherent and the phase infor- 
mation will be lost as superposition of non- coherent waves 
results in a summation of the intensities of the two waves. 
It is required that the difference between individual delays 
of electromagnetic radiation to be superpositioned is less 

15 than the wavelength of the radiation. This is a relaxed 

requirement that allows the electromagnetic radiation to be 
relatively broad-banded. For example in the visible range a 
Xe-lamp or a Hg-lamp can be used as a light source in a 
system according to the present invention with the advantage 

20 compared to a laser light source that the speckle noise is 
reduced. The requirements of the spatial coherence of the 
electromagnetic radiation depend upon the space bandwith 
product of the corresponding system and how close the 
required system performance is to the theoretically obtain- 

25 able performance of the system. 



Preferably, the electromagnetic radiation is generated by a 
coherent source of electromagnetic radiation, such as a 
laser, a maser, a phase-locked laser diode array, etc. How- 
ever a high pressure arc lamp, such as a Hg lamp, a Xe lamp, 
30 etc, may also be used and even an incandescent lamp may be 
used as a source of electromagnetic radiation in a low per- 
formance system. 



A spatial phase mask is a component that changes the phase of 
an electromagnetic wave incident upon it. The spatial phase 
35 mask may transmit or reflect the incident electromagnetic 
wave. Typically, the spatial phase mask is divided into a 
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number of resolution elements each of which modulates the 
incident electromagnetic wave by changing its phase by a 
specific predetermined value. The predetermined values are 
assigned to each resolution element in different ways depend- 
5 ing upon the technology applied in the component. For example 
in spatial light modulators, each resolution element may be 
addressed either optically or electrically. The electrical 
addressing technique resembles the addressing technique of 
solid-state memories in that each resolution element can be 

10 addressed through electronic circuitry to receive a control 
signal corresponding to the phase change to be generated by 
the addressed resolution element. The optical addressing 
technique addresses each resolution element by pointing a 
light beam on it, the intensity of the light beam correspon- 

15 ding to the phase change to be generated by the resolution 
element illuminated by the light beam. 

Spatial phase masks may be realized utilizing fixed phase 
masks, devices comprising liquid crystals and being based on 
liquid crystal display technology, dynamic mirror devices, 
20 digital micromirror arrays, defonnable mirror devices, mem- 
brane spatial light modulators, laser diode arrays (inte- 
grated light source and phase modulator), smart pixel arrays, 
etc . 

A spatial phase filter is typically a fixed phase mask, such 
25 as an optically flat glass plate coated with a dielectric 

layer at specific positions of the glass plate. However, the 
spatial phase masks mentioned in the previous section may 
also be used for spatial phase filters. 

The imaging system maps the phase modulating resolution 
3 0 elements of the spatial phase mask on the target surface of 
the synthesized intensity pattern. It may comprise a 4f-lens 
configuration (two Fourier transforming lenses utilizing 
transmission of light or one Fourier transforming lens util- 
izing reflection of light) or a single imaging lens. However, 
35 any optical imaging system providing a filtering plane for 
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the spatial phase filter may be applied in a phase contrast 
imaging system. 



In the method according to the present invention, the syn- 
thesized intensity pattern is generated by superposition of 
5 two electromagnetic waves in the image plane of the imaging 
system. The spatial phase mask changes the phase values of an 
electromagnetic wave incident upon it and the imaging system 
directs the electromagnetic wave with changed phases 
reflected from or transmitted through the spatial phase mask 

10 towards the spatial phase filter. The phase filter phase 

shifts a part of the electromagnetic radiation and the imag- 
ing system is adapted to superimpose in the image plane the 
phase shifted part of the electromagnetic radiation with the 
part of the electromagnetic radiation that is not phase 

15 shifted by the spatial phase filter. 

According to a preferred embodiment of the invention, the 
spatial phase mask is positioned at the front focal plane of 
a lens while the spatial phase filter is positioned in the 
back focal plane of the lens, whereby a first electromagnetic 

2 0 field at the phase mask is Fourier transformed by the lens 
into a second electromagnetic field at the phase filter. 
Thus, specific spatial frequencies of the first electromag- 
netic field will be transmitted through the spatial phase 
filter at specific positions of the phase filter. For 

2 5 instance, the energy of the electromagnetic radiation at zero 
frequency (DC) is transmitted through the phase filter at the 
intersecting point of the Fourier plane and the optical axis 
of the lens also denoted the zero- order diffraction region. 

It is presently preferred that the spatial phase filter is 
30 adapted to phase shift the DC-part of the electromagnetic 

radiation and to leave the remaining part of the electromag- 
netic radiation unchanged or, alternatively, to leave the DC- 
part of the electromagnetic radiation unchanged and to phase 
shift the remaining part of the electromagnetic radiation. 
35 The last alternative is preferred when the energy level of 
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the DC-part of the electromagnetic radiation is so high that 
the phase shifting part of the phase filter will be destroyed 
by it. For example in laser cutting, the DC level of the 
laser beam can be so high that a phase shifting dot posi- 
5 tioned at the intersecting point of the DC part of the laser 
beam at the phase filter would evaporate. It is also possible 
to block the electromagnetic radiation (no transmittance) in 
the zero-order diffraction region, however, the DC energy of 
the radiation is then lost. 

10 Below, an expression of the intensity of the synthesized 

intensity pattern as a function of the phasor values 4>(x,y) 
of the phase mask, when the DC-part of the electromagnetic 
radiation is phase shifted, is deduced. 

Electromagnetic radiation incident on the spatial phase mask 
15 can be described by a function A(x,y) , where A(x,y) is a 

complex number (amplitude and phase) of the incident field on 
the point (x,y) of the spatial phase mask. At the point 
(x,y) , the spatial phase mask modulates the phase of the 
incident radiation with a value <f>(x,y) so that the field 
20 after reflection by or transmission through the spatial phase 
mask may be described by the function A(x,y) * e 1 *^'*^, 
e ±<f>(x,y) b e i n g the phasor value of the point (x,y) of the 
spatial phase mask. As A(x,y) preferably is a constant value 
over the entire surface of the spatial phase mask, the term 
25 is left out of the following equations for simplicity. 

The expression of the electromagnetic radiation incident on 
the spatial phase filter may now be separated into an AC- term 
and a DC- term. If the DC- term of the field is denoted Of, the 
AC-term of the field is given by the term e iof(x ' y) As the 
3 0 spatial phase filter changes the phase of the DC-part of the 
electromagnetic radiation by 0, the intensity of the syn- 
thesized intensity pattern at the image plane of the imaging 
system is given by: 

Kx',y') = le^tx'-y'* + a(e ie -l)\ 2 
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wherein (x',y') is the coordinates of the image of the point 
(x,y) of the spatial phase mask formed by the imaging system 
in the image plane. 

It should be noted that the second term of the equation is a 
5 complex number that adds to the phasors e 1 *'*'^ of the spa- 
tial phase mask and may be interpreted as a contrast control 
parameter for the synthesized intensity pattern i(x',y'). 

According to a preferred embodiment of the invention, the 
average value of the phasors is adjusted in order to control 
10 the range of intensity levels. 

instead of phase shifting the DC-part of the electromagnetic 
radiation, it is also possible to synthesize a prescribed 
intensity pattern by phase shifting other parts of the elec- 
tromagnetic radiation by adapting the spatial phase filter to 

15 phase shift electromagnetic radiation incident upon one or 
more arbitrary regions of the phase filter and leaving the 
phase of the remaining part of the electromagnetic, radiation 
unchanged and then superimposing the two parts of the elec- 
tromagnetic radiation. The corresponding mathematics and the 

20 corresponding design procedures for the spatial phase mask 
and spatial phase filter will of course be more complicated 
than for the method described in the previous section. 

A simple example of phase shifting a part of the electromag- 
netic radiation of a spatial frequency different from the 

25 zero frequency is provided by moving the DC-part of the 

electromagnetic radiation to another spatial frequency in the 
Fourier plane (identical to the plane of the spatial phase 
filter) utilizing an optical component with an appropriate 
carrier frequency (i.e. a grating or a prism) or, preferably, 

30 encoding the function of a grating or a prism into the spa- 
tial phase mask, and adapting the spatial phase filter to 
change the phase of the electromagnetic radiation at this 
spatial frequency and to leave the phase of the remaining 
part of the electromagnetic radiation unchanged. 
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According to another preferred embodiment of the invention, 
the phase mask is not positioned in the back focal plane of 
the lens but in the Fresnel region of the lens instead. In 
this case, the electromagnetic field at the phase filter will 
5 be given by a Fresnel transformation of the electromagnetic 
field at the spatial phase mask. This further complicates the 
mathematics and the design procedures, for example the term <* 
in equation (7) has to be substituted by the value of the 
Fresnel transformation at the point (s) of phase changes of 
10 the phase filter. However, the Fresnel transformation may be 
calculated from a Fourier transformation by multiplication of 
the phasor values of the spatial phase mask by a quadratic 
phase factor followed by a Fourier transformation. 

It is an important aspect of the present invention that each 
15 intensity level of the synthesized intensity pattern for each 
resolution element may be generated by at least two different 
phasor values of a resolution element of the spatial phase 
mask. 

For example, when the spatial phase filter phase shifts the 
20 DC-part of the electromagnetic radiation, it will be shown 
later that, advantageously, the average o? of the phasors of 
the resolution elements of the phase mask should be equal to 
# and the value of the phase shift 0 should be equal to tt. In 
this case, the intensity of the synthesized image pattern at 
25 the image (x',y') of the resolution element (x,y) will be 
given by: 

I(x',y') . 2(1 - cos</>fx',y;) 

It is seen that complex conjugate phasors (values of <f> of 
opposite sign) result in identical intensity levels I(x',y'). 
30 It can be shown that for any value of the modulus of the 
average of the phasors | a| , two phasors exist that will 
generate identical intensity levels of the synthesized inten- 
sity pattern. 
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Further, if the spatial phase filter phase shifts parts of 
the electromagnetic radiation different from the DC-part , the 
phasor value that generates a specific intensity level will 
depend on the position of the resolution element in question, 
5 i.e. the phasor value and the position of the resolution 

element with that phasor value together define the intensity 
level at the image of the resolution element in the syn- 
thesized intensity pattern. Still, it is true that for each 
resolution element of the spatial phase mask, each intensity 
10 level of the synthesized intensity pattern may be represented 
by one of two different phasors of complementary phase 
values . 

This freedom of being able to select, for each intensity 
level to be generated and for each resolution element of the 

15 spatial phase mask, one of two phasors is used to control the 
phase of the Fourier transform of the phasors at specific 
spatial frequencies by selection of phasors with appropriate 
phase values to ensure two intervals of biunique functional 
dependence between phasor values and corresponding intensity 

20 values. 

This freedom of choice of phasors may be utilized to select 
phasors of neighbouring resolution elements of the spatial 
phase mask with a maximum difference between them, thereby 
generating an electromagnetic radiation emitted from the 

25 phase mask with a maximum content of high spatial frequencies 
which will generate a good separation of the DC part of the 
electromagnetic radiation from its AC part. However, any 
other strategy of selecting between two possible phasor 
values of each resolution element may be chosen to generate a 

30 desired spatial frequency content of the electromagnetic 
radiation. 

Preferably, the phase of the Fourier transform of the phasors 
at specific spatial frequencies is adjusted in order to 
control whether the relation between each phasor and the 
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corresponding intensity level is a monotonic increasing or a 
monotonic decreasing function. 

Below, a set of different methods are described that are 
provided according to the present invention for adjustment of 
5 the modulus of the Fourier transform of the phasors at speci- 
fic spatial frequencies to attain a prescribed value. If 
convenient, the methods may be combined. 

According to one of the methods, the individual phasors of 
the resolution elements of the phase mask are adjusted by a 
10 constant value until the desired value of the modulus of the 
Fourier transform of the phasors at specific spatial fre- 
quencies is attained while maintaining prescribed relative 
intensity levels between intensities of resolution elements 
of the intensity pattern, i.e. iteratively. 

15 According to another method, the individual phasors of the 
resolution elements of the phase mask are adjusted utilizing 
histogram techniques known from image processing. A histogram 
is a bar chart showing the number of resolution elements of 
the synthesized intensity pattern with a specific intensity 

20 value as a function of the intensity value. Any histogram 
technique, such as histogram equalization, adapting the 
histogram to a predetermined distribution, etc., may be used 
iteratively until the modulus of the Fourier transform of the 
phasors at specific spatial frequencies attain the prescribed 

25 value. 

According to yet another method, the phasor pattern of the 
phase mask is spatially scaled in order to adjust the modulus 
of the Fourier transform of the phasors at specific spatial 
frequencies . 

30 According to still another method, the modulus of the Fourier 
transform of the phasors at specific spatial frequencies is 
adjusted utilizing half tone coding techniques, such as 
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raster techniques, area ratio modulation, spot diameter 
modulation, etc. 

It is seen from the description above that the intensity 
levels may differ from one synthesized intensity pattern to 
the next as a consequence of the adjustments of the modulus 
of the Fourier transform of the phasors at specific spatial 
frequencies. Thus, it is preferred to control the power of 
the radiation source in dependence of the intensity range of 
the intensity pattern so that a sequence of different inten- 
sity patterns show uniform intensity levels. 

According to a preferred embodiment of the invention, the 
shape of the phase filter is adapted to match the spatial 
frequency content of the phasors of the spatial phase mask, 
e.g. to optimize the desired separation of the part of the 
electromagnetic radiation to be phase filtered from the 
remaining part of the electromagnetic radiation . 

It is within the scope of the present invention that the 
imaging system further comprises zooming means for variable 
scaling of the synthesized intensity pattern. The zooming of 
the imaging system may be dynamically controllable, e.g. in 
response to the scaling of the pattern of phasor values of 
the phase mask. 

According to the present invention, the power of the radi- 
ation source may be controllable in response to the spatial 
scaling of the pattern in the phase mask and/or the zooming 
of the focusing system. 

In order to provide a compact and integrated system according 
to the present invention, the optical function of a Fourier- 
transforming lens is encoded into the phasors of the spatial 
phase mask. The Fourier transforming lens may be refractively 
or diffract ively encoded into the phase mask. 
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Similarly, the optical function of an output lens may be 
encoded into the phase filter either refractively or 
diffract ively. 

Further, a compensation may be encoded into the phasor values 
5 of the spatial phase mask so that part of the electromagnetic 
radiation modulated by the phase mask has a substantially 
flat intensity profile in the image plane. Without this 
compensation, part of the electromagnetic radiation modulated 
by the phase mask will have a flat profile with perturbations 
10 resulting from the phase filtering superpositioned upon it. 
This may cause " ringings (oscillations) at the edges of the 
synthesized intensity pattern. 

According to another preferred embodiment of the invention, 
the source of electromagnetic radiation comprises one or more 

15 light sources of different wavelengths corresponding to three 
different colours, such as red, green and blue, for gener- 
ation of intensity patterns of arbitrary colours. Further, 
several independent systems each one illuminated by its own 
wavelength can be combined into a single mult i -wavelength 

20 system. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows a 4f common path interferometer, 

Fig. 2 shows a 2f common path interferometer. 

Fig. 3 shows a If common path interferometer, 

25 Fig. 4 shows (A) off -axis read- out of reflective SLM and (B) 
on-axis read- out of reflective SLM. 

Fig. 5 shows schematically an example of a prescribed inten- 
sity pattern in ID. 
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Fig. 6 shows schematically the resulting phase encoding 
corresponding to Fig. 5. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 



Fig. 1 shows a 4f common path interferometer (l) . A laser (2) 
5 emits a light beam which is expanded by a beam expander (3) 
into a plane light wave of uniform intensity and directs it 
towards a combination of an encrypted mask and a complex 
spatial electromagnetic radiation modulator (4) . The light 
beam is transmitted through the combination (4) and a Fourier 
10 transforming lens (5) . The combination (4) is positioned in 
the front focal plane of the lens (5) and a spatial filter 
(6) is positioned in the back focal plane of the lens (5) 
that is also the front focal plane of a lens (7) . The Fourier 
transforming lenses (5, 7) need not have identical focal 
15 lengths. Different focal lengths lead to a magnification 
ratio different from one. The filter (6) phase shifts by 0 
and optionally attenuates (by a factor B) the zero order 
diffraction part (8) of the light modulated by the 
combination (4) . Optionally, the remaining diffraction part 
2 0 of the light modulated by the combination (4) may be 

attenuated by a factor A. The synthesized intensity pattern 
is generated in the back focal plane (9) of the lens (7) and 
a dynamic focusing system (10) images the synthesized 
intensity pattern onto a focusing plane (ll) . 



25 The optical system is controlled by a computer (12) . The 

computer (12) comprises interface means for addressing each 
of the resolution elements of the complex spatial 
electromagnetic radiation modulator (4) and transmitting a 
decrypting value to the addressed resolution element. 

30 Optionally, the phase shift 0 and attenuation factors (A, B) 
of the filter (6) is adjustable and controllable by optional 
phase control means of the computer (12) which may be further 
adapted to adjust the phase shift, e.g. utilizing equation 
(18) . 
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Fig. 2 shows a 2f common path interferometer (20) . A laser 

(21) emits a light beam which is expanded by a beam expander 

(22) into a plane light wave of uniform intensity and directs 
it towards a combination of an encrypted mask and a complex 

; spatial electromagnetic radiation modulator (23) and a 
polarization beam splitter (24) and a quarter-wave plate 
(25) . The polarization beam splitter (24) and the quarter- 
wave plate (25) allows beam- splitting of light of a specific 
linear polarization without the power loss associated with 
conventional beam-splitters due to splitting of the beam in 
both directions of transmission through the beam-splitter. 
After transmission through the polarization beam splitter 
(24) and the quarter-wave plate (25) , the light beam is 
transmitted through a Fourier transforming lens (26) and is 
reflected from a spatial filter (27) . The combination (23) is 
positioned in the front focal plane of the lens (26) and the 
spatial filter (27) is positioned in the back focal plane of 
the lens (26) . The filter (27) phase shifts by B and 
optionally attenuates (by a factor B) the zero order 
diffraction part (28) of the light modulated by the 
combination (23). Optionally, the remaining diffraction part 
of the light modulated by the combination (23) may be 
attenuated by a factor A. The synthesized intensity pattern 
is generated in the back focal plane (29) of the lens (26) 
and a dynamic focusing system (3 0) images the synthesized 
intensity pattern onto a focusing plane (31) . As described 
for the system shown in Fig. l, the system (20) is controlled 
by a computer (32) . 

Fig. 3 shows a If common path interferometer (40) . A laser 

(41) emits a light beam which is expanded by a beam expander 

(42) into a plane light wave of uniform intensity and directs 
it towards a combination of an encrypted mask and a complex 
spatial electromagnetic radiation modulator (43) . The light 
beam is transmitted through the combination (43) and an image 
forming lens (44) . The filter (45) phase shifts by 6 and 
optionally attenuates (by a factor B) the zero order 
diffraction part of the light modulated by the combination 
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(43) - Optionally, the remaining diffraction part of the light 
modulated by the combination (43) may be attenuated by a 
factor A. The synthesized intensity pattern is generated in 
the image plane (46) of the lens (44) and a dynamic focusing 
5 system (47) images the synthesized intensity pattern onto a 
focusing plane (48) . As described for the system shown in 
Fig. 1, the system (40) is controlled by a computer (49) • 

Fig. 4 shows details of (Fig. 4A) an off -axis read -out of a 
combination of a reflective mask and a complex spatial 
10 electromagnetic radiation modulator (50) and of (Fig. 4B) an 
on-axis read-out of combination of a reflective mask and a 
complex spatial electromagnetic radiation modulator (51) with 
a beam splitter (52) . Both configurations (Fig, 4A, Fig. 4B) 
may be utilized in the systems shown in Figs. 1-3. 
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PHASE ENCODING FOR DC PHASE FILTERING 



In the following an example of encoding a spatial phase 
mask and a spatial phase filter will be given based on a 
5 system filtering in the DC-frequency range. The system 
chosen in this example is based on a 4-f lens 
configuration as shown in Fig, 1 and illuminated by 
electromagnetic radiation in the visible frequency 
domain, hereafter simply denoted as light radiation. 



Assuming that the illuminating light is monochromatic and 

has a substantially flat amplitude profile we obtain the 

following spatial amplitude distribution emitted from the 
spatial phase mask: 



where a(x, y) = exp(i^(x, y)) represent the spatially encoded 
phasor values and AxAy is the area of the input phase 
20 modulating spatial light modulator. 

It turns out to be convenient to separate a(x, y) into two 

terms describing a spatially invariant DC-value, a, and a 
spatially varying AC-contribution Aa(x, y) . The DC-value 
25 can be found as: 



10 



15 




(1) 




(2) 



Subsequently the AC-term is expressed by: 



30 



Aa(x, y) = exp(i#x, y)) - 



1 



JJexp(i#x,y))dxdy (3) 



AxAy 




The separation of a(x, y) into a spatially invariant DC- 
term and a spatially varying AC-term is an important 
point and will be used throughout the remaining part of 
5 this example, especially in the description of the 
spatial filtering procedure. 

The spatial filter utilized in this example is chosen as 
a circular phase contrast filter (different transverse 
10 shapes can also be used) centered around origo in the the 

spatial frequency domain, denoted by coordinates (f x , f y ) : 

7{f r ) = 1 + (exp(i0) - l)circ[-^i (4) 



15 where f z = ^Jf* + f* denotes radial spatial frequency and 

Af r describes the size of the circular (circ) phase 
filter. 

In the spatial frequency domain (the filtering plane) the 
20 Fourier transformation (3) of the spatially modulated 
light radiation from the spatial phase mask is present. 
The filtering operation on the Fourier transformed light 
radiation performed by the spatial phase contrast filter 
can be expressed as a simple point-by-point 
25 multiplication procedure. Subsequently the spatially 

filtered light is inverse Fourier transformed (3T 1 ) by the 
second Fourier lens (Fourier transformation and reflected 
output coordinates) and the resulting spatial amplitude 
distribution in the image plane (with coordinates (x f , y') ) 
30 can accordingly be written as: 



o(x', y') = a(x', y') + (exp(i0) - l)3^3(a(x, y)) c irc^Jj 

= a(x' f y') + a(exp(i0) - l)rect^ , (5) 
= [exp^x', y')) + a(exp(i0) - i)]rect^£- f ^ 



Within the illumination-region, (x', y') € 91' , outlined by 

rectf , -7— I , one obtains: 
V^Ax Ay ^ 

5 



|o(x', yf s 1 + fl sin (f)[R sin (f ) " s±T {h " * x '' ^) + f )] (6 > 

Requiring that |o(Xo' yo)f = 0 corresponding to complete 
darkness as the lowest intensity level in 
10 regions (x^, y' a ) e 9!£ implies: 

1 + sin^ljpl sin(|) - slifa - *> 0 + I)] = 0 (7) 

where the abbreviation ^ 0 = ^(x^, y^) has been used. 

15 

The solutions to Eq. (7) are given by: 



2 sir 



20 



The requirement o < |o| < 1 implies that: 



6 = n - 2(^ 5 - + p2/r, p = 0,±1,±2, 



leading to 



10 and the -sign is for ^-values: 



0 e P ; — 
J3 3 



The corresponding interval for - is: 



15 



(9) 



5 

where the +sign is for 0-values in the interval 



e e J| ; ^[ + p evan 2n (11) 



+ p^n (12) 



(* 5 -*.)e]f.-f[ 113) 

Inserting the expression for |a| / one obtains the simple 
intensity expression: 

20 

|o(x', y')f - 2^1 T sin^- - ^x', y') + |jj (14) 



where 



JJ exp(i^x, y))dxdy - AxAyp| exp(i^-) (15) 



AxAy 



25 The phase-only transformations imply that energy is 
conserved: 



JJ |o(x', y'fdx'dy' = jj\a(x t y)fdxdy = AxAy (16) 



AxAy AxAy 



5 A special case: 

The most convenient choice for a is: a = ^ (implying 

that 0 = 7T + P mn 2x) f so that the output intensity can be 
described as: 



10 |o(x', y')f = 2[l - cos(#x', y'))] (17) 

In this case the phase— »intensity mapping is described by 
the intervals [0; x] — > [0;4] . 

By setting a = i one obtains the following requirement to 
15 the phase function tf(x, y) : 



JJ cos(#x, y))dxdy = 

1r r , x (18) 
J J sin^x, y)pxdy = 0 

k AxAy 



20 



Inserting the expression for |o(x', y')f in Eq. (16) yields: 



2 III 1 ~ cos(#x', y'))]dx'dy' = AxAy (19) 

AxAy 

in accordance with the first of the integral expressions 
in Eq. (18) . 



25 



Encoding procedure: 

• A given intensity distribution (image) |o(x', y')f is 
desired at the output side of the optical setup . 

5 

• Pixellation of the image, that is generally 
represented in the greyscale range: [0; gmax] , provides 
the relation: 

gmax, 

\<\i* J)\ = 



10 



15 



25 



JJ |°<x', yfdx'dy' = AxAy ^ £ i>f = P ix i^y) 



The histogram for the desired image |o(i, j)| 2 is 

adjusted (adj) within the greyscale range [0; gmax] , so 
that the previous point is fulfilled: 

|o(i, jf |o(i, 



•c 



The phase values can now be calculated as: 
4(1, j) = arccos 



gmax 



J 



• As before pixellation provides the relation: 
20 2 sin Wi J)) = 0 • 



The previous point can now be fulfilled by complex 
conjugating half the input pixels having the same 
phase value in the phase histogram. 

The phase conjugate phase flipping provides a valuable 
tool (an extra degree of freedom) for manipulating the 
spatial frequency content in order to optimize the 



separation of low and high frequency terms at the 
filter plane. 

• The scheme is robust to constant phase errors across 
the input spatial phase modulator, since Eq. (14) is a 
function of the difference: ^ - j) , only. 
Furthermore, small variations in the individual pixel 
phase values do not introduce any detrimental effects 

because the average value a , is a result of a very 
large phasor sum. 

• If the desired intensity distribution is too small to 
include all energy, that is, the histogram is scaled 
to maximum and the left hand side of Eq. (16) is still 
smaller than the right hand side, then the input phase 
object can be scaled until Eq. (16) is fulfilled. In 
order to obtain a scale invariant output intensity 
level a dynamic focusing system is needed. Similarly, 
intensity invariance can be obtained by controlling 
the radiated power from the light source. 
Alternatively, one can ignore the residual background 
illumination and obtain intensity levels with a gain 
factor of 9*" (background constant equal to 1") for 
narrow generally shaped line structures (e.g. Eq. 
(6)). 

Exampl el: 

A very simple example illustrating the individual steps 
in the above procedure will be given below. To simplify 
the example it will be considered in one dimension only. 
The starting point for encoding the spatial phase mask in 
this example is based on the following parameters: 



10 



a = 0.5 

9 - " (20, 

gmax = 4 



Consider the pixellated 3-step function shown in Fig. 5 
to be synthesized in the image plane as an intensity 
distribution. From the above choices of parameters one 
obtains the simple relation between phase values in the 
spatial phase mask and the image intensity values: 

|o(i)| 2 = 2[l - cos(rfi))] (21) 

To proceed from here it necessary to calculate the 

accumulated intensity £ loOOf ^ n the image to be 

i 

synthesized. The accumulated intensity is easily 
calculated from an image histogram where the x-axis 
15 represents greylevel value and the y-axis represents the 
amount of pixels in the image at a given greylevel value. 

By use of a histogram ^ jo(i)| 2 is simply found as the 

i 

weighted sum of all greylevel values (x-axis) multiplied 
by their pixel counting (y-axis) . This describes, so to 
20 speak, the "weight" of the image. In this simple example 
histogram calculations are not needed since we only have 
3 greylevels with well-defined separations. 

The value for the accumulated intensity has to obey the 
25 equality: 

ZK-Ol 2 = ^#pix (Ajt) =#pix (Ax) = 14 (22) 



From Fig. 5 we obtain: 




2 \c(i)f = 4pixels • 0 + 4pixels • (0.5max) + 6pixels • max - Bmax 
j 

(23) 

5 So that the value for max can be estimated to be: 

7 

max = - (24) 
4 

The corresponding adjusted intensity levels, jo^)^, are 

10 therefore: 7/4, 7/8 and 0. These values can now be 
utilized to calculate the phase values of the spatial 
phase mask from the relation: 




(25) 



where from we obtain the three phase values: 1.4 5 rad. 
0.97 rad. and 0 rad. 

The last step needed in order to encode the spatial phase 
mask is that the following equality is fulfilled: 

20 

£sin(^(i)) = 0 (26) 
j 

Since we have the choice to use complex conjugate phasor 
values (two phasor s giving the same intensity level) many 
25 approaches can be taken from here. A simple approach is 
to flip every second phasor with its complex conjugate 
value as shown in Fig. 6. The final phase values used in 
the phase mask are accordingly: ± 1.45 rad. ± 0.97 rad. 
and 0 rad. 



30 




As the last step we can check whether the criteria: 

a = 1/2/ is actually fulfilled with the chosen phasor 
encoding: 

5 a = — (4 exp(iO) + 2 exp(i0.97) + 3 exp(£L45) + 2 exp(-i0.97) + 3 exp(-il.4 5)) = 1/2 
14 

(27) 





GENERAL PHASE .CORRECTION PROCEDURE INTEGRATED WITH THE 
PHASE ENCODING 

In Eq. (6) we obtained an analytic relation between the 
5 phase values in the spatial phase mask and the resulting 
intensity distribution, within the region (x', y') e 9?' : 



10 The analysis leading to the above relation was based on 

the assumption that [oj is a constant value within the 9T- 

domain. In other words, the following approximation was 
applied: 



However, for certain spatial filter parameters the 
lefthand side of this expression will not be a space 
invariant constant value throughout the whole 9?' -domain 
but will instead manifest slowly variations/oscillations* 

20 This will introduce small errors in the final 
superposition between the phase filtered DC-value and the 
direct propagated AC-signal. In order to circumvent this 
problem a technique is needed that can counteract the 
distortions by use of phase-only encoding in the 

25 components already present in the system. In what follows 
a procedure for integrating predistortion that 
counteracts the above mentioned distortions will be 
described that is purely based on modifying the phasor 
values in the spatial phase mask at the input side of the 

30 system. The method can also counteract other types of 
distortions inherent in a practical implementation of the 
system. Furthermore, the method can be applied in systems 
filtering at other spatial frequencies than DC. 




•Jj (28) 
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Procedure: 

When encoding the input phase function it is helpful to 
have a 'reverse' equation, expressing the input phase 
5 distribution as a function of an adjusted (electronic) 
image grey-level distribution, I sla , addressing the input 
spatial light modulator: 



15 



20 



25 



4J * Jg = 1 + 



4^(x', y')| sin^Jpx', y'\ sing) - sinf^x', y') - <*(x', y') + |)J 

10 (30) 
where it has been taken into account that a(x', y') is not 
considered as a constant but manifests a smooth 
oscillating behaviour within the optical image domain. 
The maximum value of I 3lm is denoted gmax . 



Now, one can derive a formula for the 'grey-level 
correction' tLI sXm (x', y') that one needs to apply in order to 
encode a phase function that compensates for the spatial 
variation of the average phase value a(x', y') : 

4T -^; y0 s i + <F i*. yi sii (f I^ (x '- ^ sin (S " sin (^ (x '' y,) - « x '' y,) + 9] 



I gmax J 



(31) 

where the second relation has been derived from the first 
by setting a - 1 / 2 and $ = n . 

By inserting the second relation in the first expression 
one gets: 



(32) 



This formula is however not directly useful because it is 
related to the histogram adjusted grey-level distribution 
5 denoted by J^, . 

One needs a formula that relates the above correction 
term to the 'original' input grey-level distribution 
J(x, y) that has not been modified by histogram 
10 adjustments. This is important since the effect of the 
grey-level corrections also have to be incorporated in 
the procedure of histogram adjustments. 

The histogram scaling gives: 

15 

j(x , y) = I SJa (x,y) (33) 

where I M and I alm§im are the maximum grey-level values 
occurring in the original and the adjusted electronic 
20 grey-level distributions respectively. 

Similarly, one can apply this relation to the intensity 
correction term AJ JJiB and obtain: 

25 J(x, y) = J(x, y) + A J(x, y) = -^«- (J 9io (x, y) + A J, Jin (x, y)) 

aim, max 

(34) 

resulting in: 

30 



In order to have ^enough dynamic range' in grey-levels 
for the correction term one can derive an inequality from 
the above relation by using the fact that J MR £ gmax : 



or 



2ja m in 



Xmax ~ gmax 



max 

^ J slm, 



£ gmax 



(36) 



2(gB.in| • gmax 



(37) 



Since the first term is the dominating term in the 

10 expression for the intensity correction it will in 

practice be sufficient just to have the much simpler 
corrections: 



~. v J(x, y) 

J*™ * 2)a ffl inJ • gmax 



(38) 



15 




PROPOSED APPLICATIONS: 

-Laser machining, marking, branding, trimming, hardening, 
scribing, labeling, welding and cutting on two- and three- 
dimensional surfaces especially by use of CO2 and Nd : YAG 

laser based systems. The main advantage is that energy is 
not absorbed in the system (thereby preventing damage of 
the optical hardware) and this nonabsorbed energy is 
instead utilized to increase the intensity level of the 
desired light distribution in the image plane. High power 
can be delivered to selected regions on a work piece 
simultanously. 

-Efficient and dynamic spot-array generators based on phase 
contrast imaging. In order to provide bias or holding beams 
for arrays of optoelectronic elements, such as bistable 
elements, photonic switches and smart pixels. 

-Generation of structured light (lossless) for machine 
vision applications. E.g. periodic and skew periodic mesh 
grid illumination that can be updated in parallel. 

-Photolithographic applications (laser 3D direct writing in 
parallel without the need for sequential scanning) . E.g. 
high power laser direct writing of waveguides in Ge-doped 
silica . 

-Spatial light intensity modulation in general by use of 
pure phase modulation (radiation focusators) 

-Laser beam shaping (dynamic) . 

-Highly efficient parallel image projection without the 
need for a laser scanning device. 




-Dynamic Infrared Scene Projection (DIRSP) 
-Exposure device for grating and mask production 
-LIDAR applications 
-Laserprinting in parallel 
-Laser show applications 
-Atmosphere research 
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CLAIMS 



1. A method of decryption of an encrypted image having a non- 
encrypted image intensity pattern I(x',y') and 

encoded into a mask having a plurality of mask resolution 
5 elements (x^n,y m ) with an encoded phase value <t>(Xm,Y m ) and an 
encoded amplitude value a(x m# y m ) / and 

encrypted by addition of an encrypted phase value 4>c (Xm#y m ) 
to the encoded phase values ^fx^y^ and by multiplication of 
an encrypted amplitude value a c (x m ,y m ) with the encoded 
10 amplitude value a(x^ # y m ) , 

each mask resolution element (x^y^ modulating the phase and 
the amplitude of electromagnetic radiation incident upon it 



15 the method comprising the steps of 

radiating electromagnetic radiation towards the mask, 

inserting into the path of the electromagnetic radiation a 
complex spatial electromagnetic radiation modulator 
comprising modulator resolution elements (x d ,y d ) , each 

20 modulator resolution element (x d ,y d ) modulating the phase and 
the amplitude of electromagnetic radiation incident upon it 
with a predetermined complex value a d (x d ,y d ) e i*d<*a,yd) ^ the 
decrypting phase value <£ d (x d ,y d ) and the decrypting amplitude 
value a d (x d ,y d ) f respectively, of a modulator resolution 

25 element (x d ,y d ) being substantially equal to -<^c(x m/ y m ) and 
a c _1 {Xj,,, y m ) , respectively, of a corresponding mask resolution 
element (x^y^) , and 

imaging the mask and the electromagnetic radiation modulator 
onto the image having the image intensity pattern I(x',y'). 



with the complex value a (x^y^ a c (x^, y m ) e 
and 



i# (xm, ym) +i#c (xm, ym) 
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2. A method according to claim 1, wherein the step of imaging 
comprises imaging with a common path interferometer. 

3. A method according to claim 1 or 2, wherein the step of 
imaging comprises phase contrast imaging. 

5 4. A method according to claim 2 or 3, further comprising the 
steps of 

Fourier or Fresnel transforming electromagnetic radiation 
modulated by the mask and the complex spatial electromagnetic 
radiation modulator, 

10 filtering the Fourier or Fresnel transformed electromagnetic 
radiation by 

in a region of spatial frequencies comprising DC in the 
Fourier or Fresnel plane, 

phase shifting with a predetermined phase shift value 
15 6 the modulated electromagnetic radiation in relation 

to the remaining part of the electromagnetic 
radiation , and 

multiplying the amplitude of the modulated 
electromagnetic radiation with a constant B, and 

20 in a region of remaining spatial frequencies in the 

Fourier or Fresnel plane, multiplying the amplitude of 
the modulated electromagnetic radiation with a constant 
A, 

forming the intensity pattern by Fourier or Fresnel 
25 transforming, respectively, the phase shifted Fourier or 
Fresnel transformed modulated electromagnetic radiation, 
whereby each resolution element (Xm#y m ) of the mask is imaged 
on a corresponding resolution element (x',y') of the image, 
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the filtering parameters A, B, 0 substantially fulfilling the 
equation 

I(x',y') = A 2 |a(x',y')e i * <x ''y , > + BA^ate 16 -!) | 2 
of being the average of the complex phasors a (x^ y m ) e i * ^-V") , 

5. A method according to claim 4, wherein the filtering 
parameters A and B substantially fulfil that A=l and B«l, and 
for each (x mf y m ) of the mask: a{x mf y m )=l. 

6. A method according to claim 5, wherein the phase shift 
value 0 substantially fulfils the equation 



2 | Bin.* | 

7. A method according to claim 6, wherein the phase shift 0 
is substantially equal to 7r. 

8. A method according to any of claims 3-7, further 
comprising the steps of 

moving the DC- part of the electromagnetic radiation to a 
second part of the Fourier or Fresnel plane, and 

phase shifting the Fourier or Fresnel transformed modulated 
electromagnetic radiation at the second part of the Fourier 
or Fresnel plane by 0 in relation to the remaining part of 
the electromagnetic radiation. 

9. A method according to claim 8, wherein the step of moving 
the DC-part of the electromagnetic radiation comprises 
utilization of an optical component, such as a grating, a 
prism, etc, with an appropriate carrier frequency. 
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10. A method according to any of claims 3-9, further 
comprising the step of phase shifting at selected spatial 
frequencies constituting a region that is shaped to match the 
spatial frequency content of the phasors e i * (xm ' yin> . 

5 11. A method according to any of claims 3-10, wherein the 
step of filtering comprises utilization of a spatial light 
modulator . 

12. A method according to any of claims 3-11, further 
comprising the step of encoding the optical function of an 

10 output lens into the filter. 

13. A method according to any of claims 3-12, wherein the 
step of radiating electromagnetic radiation comprises 
radiation of electromagnetic radiation of different 
wavelengths corresponding to three different colours, such as 

15 red, green and blue, for generation of intensity patterns of 
arbitrary colours. 

14. A method of encryption of an image having an intensity 
pattern I(x',y') to be decrypted according to any of claims 
1-15, comprising the steps of 

20 pixellating the intensity pattern I(x',y') in accordance with 
the disposition of resolution elements (x^y^ of a mask, 

encoding the mask with an encoded phase value 0(x m ,y m ) and an 
encoded amplitude value afx^y^, and 

encrypting by addition of an encrypted phase value ^cfx^y^ 
25 to the encoded phase values (Mx^y^ and by multiplication of 
an encrypted amplitude value a^x^y^ with the encoded 
amplitude value a(x m ,y m ), 

each mask resolution element < 3^ , y m ) modulating the phase and 
the amplitude of electromagnetic radiation incident upon it 
30 with the complex value a (x m , y m ) a c (x m , y m ) e 1 * ( ^' ym) +i * c (xm 'y in) . 
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15. A method according to claim 14, further comprising the 
step of calculating the complex phasor values 
a < x m'ym) ei * (xin ' yin) of the mask in accordance with 

I(x',y') = A 2 laU'.yMe 1 ***''*'* + BA' 1 ^ (e ie - 1 ) | 2 

for selected phase shift values 6, a being the average of the 
complex phasors a (x^, y m ) e 1 * ^'i™) , 

selecting, for each resolution element, one of two phasor 
values which represent a particular grey level . 

16. A method according to claim 15, wherein the filtering 
parameters A and B substantially fulfil that A=l and B-l, and 
for each (x m ,y m ) of the mask: a(x m ,y m )=l. 

17. A method according to claim 16, wherein the step of 
calculating the phasor values comprises 

setting the synthesized intensity of at least one resolution 
element (x 0 ', y 0 ') of the intensity pattern to zero, and 

calculating the phasor values e d * (x ' y) of the mask in 
accordance with 

i*i = — i—ff- 

2|sin|| 

Hx'.y') = 2 [lTsin^-^U'fy') + -§)] 
for selected phase shift values 6, 4% being the phase of c*. 

18. A method according to claim 17, further comprising the 
step of selecting the phase shift 6 = ir, selecting |a| = 
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and calculating the phasor values e 1 ***'* 0 of the phase mask 
in accordance with 



19. A method according to any of claims 14-18, further 
comprising the step of encoding the function of an optical 
component, such as a grating, a prism, etc, with an 
appropriate carrier frequency, into the mask. 

20. A method according to any of claims 14-19, further 
comprising the step of adjusting the modulus of the Fourier 
transform of the complex phasors a (x^ y m ) e**'™'^ at 
specific spatial frequencies in order to control the range of 
intensity levels of the synthesized intensity pattern. 

21. A method according to claim 20, wherein the step of 
adjusting the modulus of the Fourier transform of the complex 
phasors a(x m ,y m )e i * (xin 'y iri) at specific spatial frequencies 
comprises at least one of the following measures: 

a) adjusting the individual complex phasors 
^Xm^m)^^'^ of the resolution elements of the mask 
maintaining prescribed relative intensity levels between 
intensities of resolution elements of the intensity 
pattern, 

b) adjusting the individual complex phasors 
a ( x m^ym)e x * (xin ' yin) of the resolution elements of the mask 
by histogram techniques, 

c) spatially scaling the complex phasor a (x m ,y m ) e 1 *^'^ 
pattern of the mask, and 



J(x' 



,y') = 2[l-cos(4>(x / ,y / ))] 




mask 
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d) utilizing half tone coding techniques. 

22. A method according to any of claims 14-21, wherein each 
complex phasor a <x^,y m ) e^^Y™) of the mask is selected from 
a set of two determined phasors with complementary complex 

5 phasor values a(x B ,y B )e i * 1,JBB ^» and a (x^ yj e 1 * 2 in such 

a way that a specific spatial frequency distribution of the 
intensity of the electromagnetic radiation in the Fourier or 
Fresnel plane is attained. 

23. A method according to claim 22, wherein the phase <£(x,y) 
10 of complex phasors a (x^ y m ) e i * (3an *y m) of adjacent resolution 

elements alternates between the two possible complementary 
complex phasor values a (x m ,y m ) e 1 * 1 and 
a(x m ,y m )e i * 2(3OT ^ ) . 

24. A method according to claim 22 or 23 # wherein the complex 
15 phasors a (x^yj e 1 * 1 <».Y»> and a (x^, y m ) e 1 * 2 are complex 

conjugated. 

25. A method according to any of claims 14-24, further 
comprising the step of encoding the optical function of a 
Fourier- transforming lens into the complex phasors 

20 a(x m ,y m )e i * (xin ' yin> of the mask. 

26. A decryption system for decrypting an encrypted image 
having a non-encrypted image intensity pattern I (x f ,y') that 
has been 

encoded into a mask having a plurality of mask resolution 
25 elements (x m ,y m ) with an encoded phase value *(x m ,y m ) and an 
encoded amplitude value a(x m# y m ), and 

encrypted by addition of an encrypted phase value ^ c (x m ,y m ) 
to the encoded phase values <t>(Xm,Y m ) and by multiplication of 
an encrypted amplitude value a c (x m# y nv ) with the encoded 
30 amplitude value a(x m ,y m ) , 
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each mask: resolution element (x m ,y m ) modulating the phase and 
the amplitude of electromagnetic radiation incident upon it 
with the complex value a t^, y m ) a c (x,,,, y m ) e 1 * (M ' ^ +i * c (3an ' V 1 " 1 , 

the system comprising 

5 a source of electromagnetic radiation for emission of 
electromagnetic radiation, 

positioned in the path of the electromagnetic radiation a 
complex spatial electromagnetic radiation modulator 
comprising modulator resolution elements (x d ,y d ) , each 

10 modulator resolution element (x d ,y d ) modulating the phase and 
the amplitude of electromagnetic radiation incident upon it 
with a predetermined complex value a d (x d/ y d ) e i * d(xd ' yd) , the 
decrypting phase value # d (x d ,y d ) and the decrypting amplitude 
value a d (x d ,y d ), respectively, of a modulator resolution 

15 element (x d ,y d ) being substantially equal to -^c^^'Ym^ and 
a c " 1 < x m'y , m^' respectively, of a corresponding mask resolution 
element (x m ,y m ) , and 

an imaging system for imaging the mask and the 
electromagnetic radiation modulator onto the image having the 
20 image intensity pattern I (x' ,y'). 

27. A system according to claim 26, wherein the imaging 
system comprises a common path interferometer. 

28, A system according to claim 26 or 27, wherein the imaging 
system comprises a phase contrast imaging system, 

25 29. A system according to claim 27 or 28, further comprising 

means for Fourier or Fresnel transforming the electromagnetic 
radiation modulated by the mask and the complex spatial 
electromagnetic radiation modulator and being positioned on a 
propagation axis of the modulated radiation, 
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a spatial filter for filtering the Fourier or Fresnel 
transformed electromagnetic radiation by 

in a region of spatial frequencies comprising DC in the 
Fourier or Fresnel planie, 

5 phase shifting with a predetermined phase shift value 

0 the modulated electromagnetic radiation in relation 
to the remaining part of the electromagnetic 
radiation, and 

multiplying the amplitude of the modulated 
10 electromagnetic radiation with a constant B, and 

in a region of remaining spatial frequencies in the 
Fourier or Fresnel plane, multiplying the amplitude of 
the modulated electromagnetic radiation with a constant 
A, 

15 means for forming the intensity pattern by Fourier or Fresnel 
transforming, respectively, the phase shifted Fourier or 
Fresnel transformed modulated electromagnetic radiation, 
whereby each resolution element (x^y^ of the mask is imaged 
on a corresponding resolution element (x',y') of the image, 

20 the filtering parameters A, B, 9 substantially fulfilling the 
equation 

I(x',y') = A 2 |a(x' f y')e i * (x ''y' ) + BA' 1 ^ (e ie - 1) | 2 

for selected phase shift values 6, <x being the average of the 
complex phasors a (x m , y m ) e i*(«n#Yro) # 

25 30. A system according to claim 29, wherein the filtering 

parameters A and B substantially fulfil that A=l and B=l, and 
for each (x^y^ of the mask: a(x m ,y m )=l. 
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31. A system according to claim 30, wherein the phase shift 
value 0 substantially fulfils the equation 



I*' = —77- 

2|siniL 
1 2 



32. A system according to claim 31, wherein the phase shift 6 
is substantially equal to 7r. 

5 33. A system according to any of claims 28-31, further 
comprising 

means for moving the region of spatial frequencies comprising 
DC to a second part of the Fourier or Fresnel plane, and 
wherein 



10 the spatial filter is adapted to phase shift the transformed 
modulated electromagnetic radiation at the second part of the 
Fourier or Fresnel plane by 9 in relation to the remaining 
part of the electromagnetic radiation. 



34. A system according to claim 33, wherein the means for 
15 moving the region of spatial frequencies comprising DC to a 

second part of the Fourier or Fresnel plane comprises an 
optical component, such as a grating, a prism, etc, with an 
appropriate carrier frequency. 

35. A system according to any of claims 29-34, wherein the 
20 spatial filter comprises a spatial light modulator. 

36. A system according to any of claims 29-35, wherein the 
spatial filter is adapted to perform the optical function of 
an output lens by appropriate encoding of the spatial filter. 

37. A system according to any of claims 26-36, wherein the 
25 source of electromagnetic radiation is adapted to radiate 

electromagnetic radiation of different wavelengths 
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corresponding to three different colours, such as red, green 
and blue, for generation of intensity patterns of arbitrary 
colours . 



38. A system according to any of claims 29-37, further 

5 comprising a first and a second Fourier transforming lens, 
the mask being positioned in the front focal plane of the 
first lens, the spatial filter being positioned at the back 
focal plane of the first lens, and the second lens being 
positioned so that its front focal plane is positioned at the 
10 position of the back focal plane of the first lens. 

39. A system according to any of claims 29-38, further 
comprising one Fourier transforming lens, the spatial filter 
being positioned at the back focal plane of the lens. 

40. A system according to any of claims 29-39, further 
15 comprising one imaging lens, the spatial filter being 

positioned in the back focal plane of the lens. 

41. A system according to any of claims 29-40, further 
comprising a polarizing beam splitter and a quarter wave 
plate and/or a phase filter reflecting electromagnetic 

20 radiation incident upon it. 

42. A system according to any of claims 29-41, wherein the 
spatial filter changes the phase of the radiation in the 
region of spatial frequencies comprising DC and leaves the 
phase of the remaining part of the radiation unchanged. 

25 43. A system according to any of claims 29-41, wherein the 
spatial filter do not change the phase of the radiation in 
the region of spatial frequencies comprising DC and changes 
the phase of the remaining part of the radiation. 
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44. A system according to any of claims 29-41, wherein the 
spatial filter blocks the radiation at the region of spatial 
frequencies comprising DC and leaves the remaining part of 
the radiation unchanged. 

45. A system according to any of claims 29-44, wherein the 
source of electromagnetic radiation is a Laser. 
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